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Self-consistent solutions to the nonlinear spinor field equations in General Relativity has been 
studied for the case of Bianchi type-I (B-I) space-time. It has been shown that, for some special type 
of nonliearity the model provides regular solution, but this singularity-free solutions are attained at 
the cost of broken dominant energy condition in Hawking-Penrose theorem. It has also been shown 
that the introduction of A-term in the Lagrangian generates oscillations of the B-I model, which 
is not the case in absence of A term. Moreover, for the linear spinor field, the A term provides 
oscillatory solutions, those are regular everywhere, without violating dominant energy condition. 
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1. INTRODUCTION 



Nonlinear phenomena have been one of the most popular topics during last years. Nevertheless, it must be admitted 
that nonlinear classical fields have not received general consideration. This is probably due to the mathematical 
^^ difficulties which arise because of the nonrenormalizability of the Fermi and other nonlinear couplings B . Nonlinear 
K^ [ selfcouplings of the spinor fields may arise as a consequence of the geometrical structure of the space-time and, more 
J^ I precisely, because of the existence of torsion. As soon as 1938, Ivanenko [0-0 showed that a relativistic theory imposes 
,__! ' in some cases a fourth order selfcoupling. In 1950 Weyl B proved that, if the affine and the metric properties of 
f^ [ the space-time are taken as independent, the spinor field obeys either a linear equation in a space with torsion or a 
[■ — , nonlinear one in a Reimannian space. As the selfaction is of spin-spin type, it allows the assignment of a dynamical 
f^ role to the spin and offers a clue about the origin of the nonlinearities. This question was further clarifed in some 
I important papers by Utiyama, Kibble and Sciama p-@| In the simplest scheme the selfaction is of pseudovector type, 
(•^ , but it can be shown that one can also get a scalar coupling ||]. An exelent review of the problem may be found 
|1C| ]. Nonlinear quantum Dirac fields were used by Heisenberg |11 1^ in his ambitious unified theory of elementary 






particles. They are presently the object of renewed interest since the widely known paper by Gross and Neveu \li 
X The quantum field theory in curved space-time has been a matter of great interest in recent years because oTlts 
?-H applications to cosmology and astrophysics. The evidence of existence of strong gravitational fields in our Universe 
^^' led to the study of the quantum effects of material fields in external classical gravitational field. After the appearance 
of Parker's paper on scalar fields |Q and spin-i fields, ||l5[ several authors have studied this subject. The present 
cosmology is based largely on Friedmann's solutions of the Einstein equations which describe the completely uniform 
and isotropic universe ('closed' and 'open' models, i.e., bounded or unbounded universe). The main feature of these 
solutions is their non-stationarity. The idea of expanding universe, following from this property, is confirmed by the 
astronomical observations and it now safe to assume that the isotropic model provides, in its general features, an 
adequate description of the present state of the universe. Although the Universe seems homogenous and isotropic at 
present, it does not necessarily mean that it is also suitable for description of the early stages of the development of 
the universe and there are no observational data guarantying the isotropy in the era prior to the recombination. In 
fact, there are theoretical arguments that sustain the existence of an anisotropic phase that approaches an isotropic 
one |l6||. I nterest in studying Klein-Gordon and Dirac equations in anisotropic models has increased since Hu and 
Parkerll^ have shown that the creation of scalar particles in anisotropic backgrounds can dissipate the anisotropy 
as the Universe expands. 

A Bianchi type-I (B-I) Universe, being the straightforward generalization of the flat Robertson- Walker (RW) Uni- 
verse, is one of the simplest models of an anisotropic Universe that describes a homogenous and spatially flat Universe. 
Unlike the RW Universe which has the same scale factor for each of the three spatial directions, a B-I Universe has a 
different scale factor in each direction, thereby introducing an anisotropy to the system. It moreover has the agreeable 
property that near the singularity it behaves like a Kasner Universe, even in the presence of matter, and consequently 
falls within the general analysis of the singularity given by Belinskii et al Ig^l ■ Also in a Universe flUed with matter 
for p = ^e, ^ < 1, it has been shown that any initial anisotropy in a B-I Universe quickly dies away and a B-I 
Universe eventually evolves into a FRW Universe M]. Since the present-day Universe is surprisingly isotropic, this 
feature of the B-I Universe makes it a prime candidate for studying the possible effects of an anisotropy in the early 
Universe on present-day observations. In light of the importance of mentioned above, several authors have studied 
B-I universe from different aspects. 

In p^] Chimento and MoUerach studied the Dirac equations in B-I universe and obtained their classical solutions. 
They also claimed that for each value of the momentum only two independent solutions exist and showed that it is 
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not possible to obtain the solutions from those of a FRW universe only by perturbation. One of the solutions obtained 
would describe a particle with a given helicity, while the other one would represent antiparticles with the opposite 
helicity. This fact posed a very interesting problem, spin-1/2 particles cannot live in a B-I, at least if they keep their 
well-known properties of flat space-time. This problem was handled by Castagnino et al [^ where they showed that 
if the Dirac equation is separable, the number of independent solution is four, contrary to the claim made in ||22| . 

In a number of papers p4|-|2q| several authors studied the behavior of gravitational waves (GW's) in a B-I Universe. 
In P5[ the evolution equations for small perturbations in the metric, energy density, and material velocity were 
derived for an anisotropic viscous B-I universe.lt has been shown that the results were independent of the equation of 
state of the cosmic fluid and its viscosity. They also showed that the GWs need not necessarily be transversal in an 
anisotropically expanding B-I universe and the longitudinal components of the gravitational waves have no physical 
significance. In pq] Cho and Speliotopoulos studied the propagation of classical gravitational waves in B-I universe. 
They found that CJWs in B-I universe are not equivalent to two minimally coupled massless scalar fields as in FRW 
universe. Because of its tensorial nature, the GW is much more sensitive to the anisotropy in space-time than the 
scalar field is and it gains an effective mass term. Moreover, they found a coupling between the two polarization 
states of the GW which is not present in a FRW universe. 

Nonlinear spinor field (NLSF) in external FRW cosmological gravitational field was first studied by G.N. Shikin in 
1991 p^. The main purpose to introduce a nonlinear term in the spinor field Lagrangian is to study the possibility 
of elimination of initial singularity. Following [ p7| , we analyzed the nonlinear spinor field equations in an external B-I 
universe p8{ . In that paper we consider the nonlinear term in the spinor field Lagrangian as an arbitrary function 
of all possible invariants generated from spinor bilinear forms. There we also studied the possibility of elimination of 
initial singularity especially for the Kasner Universe. For few years we studied the behavior of self-consistent NLSF 
in a B-I Universe p^,^ both in presence of perfect fluid and without it that was followed by the Refs., [3l|-p3[ where 
we studied the self-consistent system of interacting spinor and scalar fields. Recently, we study |^,^ the role of the 
cosmological constant (A) in the Lagrangian which together with Newton's gravitational constant (G) is considered 
as the fundamental constants in Einstein's theory of gravity |36[| . 

2. REVIEW OF B-I COSMOLOGY 

A diagonal Bianchi type-I space-time (hereafter B-I) is a spatially homogeneous space-time which admits an abelian 
group Ga, acting on spacelike hypersurfaces, generated by the spacelike Killing vectors Ci = 9i, ^2 = ^2, ^3 = ^3. In 
synchronous coordinates the metric is |37|j3^ : 

3 

ds^ = dt^ ~^ai{t)dxi. (2.1) 

1=1 

If the three scale factors are equal (i.e., ai — a2 = 03), ( ^.l[ ) describes an isotropic and spatially flat Friedmann- 
Robertson- Walker (FRW) universe. The B-I Universe has a different scale factor in each direction, thereby introducing 
an anisotropy to the system. Thus, a Bianchi type-I (B-I) universe, being the straightforward generalization of the flat 
Friedmann-Robertson- Walker (FRW) universe, is one of the simplest models of an anisotropic universe that describes 
a homogenous and spatially flat universe. When the two of the metric functions are equal (e.g. 02 =03) the B-I 
space-time is reduce to the important class of plane symmetric space-time (a special class of the Locally Rotational 
Symmetric space-times |lq,u9|) which admit a G4 group of isometrics acting multiply transitively on the spacelike 
hypersurfaces of homogeneity generated by the vectors ^j^, ^21 €3 ^.nd ^4 = x^d^ — x^d^- The B-I has the agreeable 
property that near the singularity it behaves like a Kasner universe, given by 

ai(t) = a\t^\ a-2.{t) = a°i''% a-i{t) = a°i2'^ (2.2) 

with pj being the parameters of the B-I space-time which measure the relative anisotropy between any two asymmetry 
axis and satisfy the constraints 

Pi+P2+P3 = l (2.3a) 

pI+pI+pI = 1- (2.3b) 

Thus out of three parameters, only one is arbitrary. One particular choice of parametrization is 

Pi - 2 r^^ 1 (2-4a) 

p'^ +P+1 



PJP + 1) 
p'^ + p+1 



p, = J^LJ-±>- (2.4b) 



^3 = ^^^. (2.4c) 
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The condition < p < 1 on p then yields the condition —1/3 < pi < 0, < P2 % 2/3, 2/3 < ps < 1. Another 



particular parametrization can be given using an angle on the unit circle, since (2.3) describes the intersection of a 
sphere with aplane in the parameter space (pi,p2,P3): 

pi = -(1 + cosi9 + TSsin??), (2.5a) 

P2 = ^(l + cos?9-\/3sim?), (2.5b) 

Pi. = -{l-2cosd). (2.5c) 

Although 1? ranges over the unit circle, the labeling of each pj is quite arbitraey. Thus the unit circle can be devided 
into six equal parts each of which span 60°, and the choice oi pj is unique within each section separately. For d — Q, 
Pi = P2 = 2/3 and ps = —1/3 while for ■& = 7r/3 pi = 1 and P2 = P3 = 0. 

Let us now go back to B-I metric. The non-trivial Christoffel symbols for (^) are 

r% = a.a,, r^, = r^o = ^, (2.6) 

while the components of non-trivial Ricci tensor read 

3 .. 

i—\ ■' 

The Ricci scalar for the B-I universe has the form 



i?oo--V-, i?,,= p-H^(^ + ^)la^ z,j,fc=l,2,3, i+3i^k. (2.7) 

"^ fli La,; ai Va, afc/J 



i?= -21 — -I- — -I- — -f — — -I-— — + — ). (2.8) 

\a\ a-i as ai 02 02 as 03 ai ) 

Sometimes it proves convenient to introduce a new time parameter t] by 

t 

r;= f a-\t)di, (2.9) 

where we define 

[ait)]^ = C{t) = (aia2a3)'/' - {CiCiC^f^ (2.10) 

with C ~ i = of. Note that in the isotropic limit, i.e., ai = ai = a^ r/ reduces to conformal time. Further, defining 



c; „ 1 V- , C" 
d 



c' ^^Je^» = ^' Q^^Y.('^.-d,r (2.11) 

i— 1 i<j 

where prime (/) denotes differentiation with respect to 77, we get the following nonzero Christoffel symbols for the 
metric ( |2.lD 

F'' --D F'' - - ' ' r^ -F* -id- ('2 12) 



The nonzero components of the Ricci tensor now read 

R^^ = Id' + 6Q, R,, = -^{d[ + d,D) (2.13) 

and the Ricci scalar 



.^^-\d' + QQ, R,, = -^('^' 



R = C-^[3D' + -D^ + 6Q]. (2.14) 



Note the in the sections to follow, we work with the usual time t. 
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3. FUNDAMENTAL EQUATIONS AND GENERAL SOLUTIONS 

The action of the nonUnear spinor, scalar and gravitational fields can be written as 

Sig;i^,i>,ip)^ J LV^dn (3.f) 

with 

L = Lg + Lsp + Ln,. (3.2) 

Here Lg corresponds to the gravitational field 

R + 2A 



2k 



(3.3) 



where R is the scalar curvature, k = 8ttG with G being the Einstein's gravitational constant and A is the cosmological 
constant. The spinor field Lagrangian Lgp is given by 






V'T^'V^t/- - v^Vt'^V' 



rmpip + Lf^, (3.4) 



where the nonlinear term i^ describes the self-interaction of a spinor field and can be presented as some arbitrary 
functions of invariants generated from the real bilinear forms of a spinor field. Since '0 and ip* (complex conjugate of 
ip) has 4 component each, one can construct 4-4=16 independent bilinear combinations. They are 

(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 
(3.5e) 



5 = ^V, 


(scalar). 


P = iipj ip^ 


(pseudoscalar). 


v^ = (ipj'^ip), 


(vector) 


A^ = {^-/^-ff'^j 


), (pseudovector). 


r^'' = (V5a^» 


(antisymmetric tensor). 


'^]. Invariants, corresponding to the bilnear foi 


/ = 5^ 




j = p', 




Iv ^V^V^ = 


{^^''i^)g^,{{pYi^), 


I A = A^ A^ = 


-- (VV7^V')ffA..(^7'7'^V'), 



(3.6a) 
(3.6b) 
(3.6c) 
(3.6d) 

It = T^. T^"^ = (V^a'^'V) g^,a.g,fi[^G'^^i>). (3.6e) 

According to the Pauli-Fierz theorem, ||40|| among the five invariants only / and J are independent as all other can 
be expressed by them: I^ — —Ia = I + J and It = I— J. Therefore we choose the nonlinear term F to be the function 
of / and J only, i.e., Ln = F{I i J)i thus claiming that it describes the nonlinearity in the most general of its form. 
Lin is the Lagr angi an of perfect fluid. 

Variation of (3.1) with respect to spinor field ip (ip) gives nonlinear spinor field equations 

ij^'S/^^ip-m^P + ViP + gi-f^ij^O, (3.7a) 

iVpV;7^ + m^/j -V^!- gii^j^ = 0, (3.7b) 

where we denote 

dF dF 

Varying (pTl) with respect to metric tensor g^y one finds the Einstein's field equation 

R>i'\5^yR^->^Tl^ + AK (3.8) 

where i?^ is the Ricci tensor; R — g^'^R^i, is the Ricci scalar; and Tj^ is the energy-momentum tensor of the material 
field given by 
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/J sp /^ ^ in /i ■ 



Here T^'^ is the energy-momentuni tensor of the spinor field 



^sp''^ = T5'"' ( ini^uip + V'Ti^v^.v - v^jVTi^V' - ^yipi,,^ 



5;ls 



where Lsp with respect to (3.7) takes the form 



L, 



-{VS + gP)+F{I,J). 
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(3.9) 

(3.10) 
(3.11) 



'^il(m) ^^ ^^^ energy- momentum tensor of a perfect fluid. For a Universe filled with perfect fluid, in the concomitant 
system of reference (u*^ = 1, u* = 0, i = 1, 2, 3) we have 



r/^(m) = (P + £)u^,u" - S^^p = (e, ~p, -p, -p), 



(3.12) 



where energy £ is related to the pressure p by the equation of state p — Qe. The general solution has been derived 
by Jacobs |£l|. Here C varies between the interval < ^ < 1, whereas C = describes the dust Universe, C, — \ 
presents radiation Universe, -^ < C < 1 ascribes hard Universe and C = 1 corresponds to the stiff matter. 

In (3.7) and (3.9) V^ denotes the covariant differentiation; its explicit form depends on the quantity it acts on. 
This covariant differentiation has the standard properties 



V^7. = 0, 



(3.13a) 
(3.13b) 
(3.13c) 



where the symbol * means Hermitian adjoint (the transpose of the complex conjugate). The explicit form of the 
covariant derivative of spinor is pd| , [4^ 



dip 



(3.14a) 
(3.14b) 



where r^(x) are spinor affine connection matrices. 7 matrices in the above equations obey the following algebra 



7^7"^ + 7''7^ = 2g'"' 
and are connected with the flat space-time Dirac matrices 7 in the following way 

gtj.,.{x) = e''{x)el{x)i]ab, l^,{x) = e°(x)7Q, 



(3.15) 



(3.16) 



where rjab — diag(l, —1, —1, —1) and e" is a set of tetrad 4-vectors. The spinor afhne connection matrices r^(x) are 
uniquely determined up to an additive multiple of the unit matrix by the equation 



^^71^ = -g^ - rg^7p - r^7^ + 7^r^ = o. 



r,(x) = -5p.(x)(VK-r^j7V- 



with the solution 



Let us now write the 7's and F^'s explicitly for the B-I metric ( p.l[ ) that we rewrite in the form 

ds^ ^ dt^ - a^{t)dx^ - b^{t)dy^ - c^{t)dz^. 
For the metric ( 3.19 ) from ( |3.16| ) one finds 

70 = 70, 71 = a{t)^i, 72 = 6(i)72, 73 = 2(^)73, 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



7" =7", 7'=7V«W, l^^flKt), I'^f/cit). 
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For the affine spinor connections from (3.1q) we find 



ro = o, Ti = id(t)7'7°, r2 = ife(i)7V, r3 = ic(i)f7°- 



Flat space-time matrices 7 we will choose in the form, given in 

10 
10 
0-10 
0-1 



-i 

i 

i 

-i 



f = 



/ 1 
10 
0-100 

V-1 0. 

10 
0-1 
-10 
10 



Defining 7^ as follows, 



7' = -^i^M-p7''7''7''7^ 
7^ = -«\/^7°7S^7^ = 



E^uap — V^g^l-ivapj £0123 — 1) 



.i^O-1-2-3 ^ -5^ 



we obtain 



f = 



/ 0-10 
0-1 
-10 
\ -1 



For the space-time (3.19) Einstein equations (3.8) now read 



(3.21) 



b c be 
b e be 


^kTI 


-A, 


e a ea 
e a ea 


= nTi 


-A, 


a b ab 

- + I + -I: 
a b ab 


= KTi 


-A, 


ab be ea 

-7 + 7- + -- 
ab be ea 


= -r„" 


-A, 



where point means differentiation with respect to t. 



We will study the space-independent solutions to the spinor field equations (3.7) so that ip — tp{t). Setting 

T — abe = \/— g 



we rewrite the spinor field equation (3.7a) as 



Z7 



d 



dt 2t 



+ --?/;" "11/' + 'Dil} + Gij^ip = 0. 



Setting Vj{t) — ^/Tipj{t), j = 1,2,3,4, from (3.24) one deduces the following system of equations 



Vi+i{m-V)Vi-gV3=0, 
V2+i{m-V)V2-gVi^0, 
V3^iim-V)V3+gVi=0, 

Using the solutions obtained one can write the components of spinor current: 



(3.22a) 
(3.22b) 

(3.22c) 

(3.22d) 

(3.23) 
(3.24) 



(3.25a) 
(3.25b) 
(3.25c) 
(3.25d) 
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jA* = ^^M^_ (3,26) 

Taking into account that ^ = "0^7°, where -0^ = {^l, -02, "03: i^l) ^-nd V'j = ^j/V^j J — I7 2, 3, 4 for the components 
of spin current we write 

f = - [V^Vi + V;Vi + ^3*1^3 + ^4*^4] , (3.27a) 

3^ = — \V{V4 + V^Vj. + ViV2 + V^Vi] , (3.27b) 

ar 

f = ^ [^i^i - V2V^ + ViV2 - VlVi] , (3.27c) 

/ = ^ [^1*^3 - V^Vi + V^Vi - V:V2] . (3.27d) 

The component j° defines the charge density of spinor field that has the foUowing chronometric-invariant form 

e=(jo-J°)'/^ (3.28) 

The total charge of spinor field is defined as 



Q= I gy^-^gdxdydz (3.29) 

Let us consider the spin tensor [Q 

S^""'' = ^1^(7^^" + a^^'Y}^- (3.30) 



We write the components 5"*^'° (z,fc = 1,2,3), defining the spatial density of spin vector explicitly. From ( 3.3C ) we 
have 

S-^'O = ^^P{i°cf'^ + a'^-f°}tP = i07V'^ V (3.31) 

that defines the projection of spin vector on k axis. Here i, j, k takes the value 1, 2, 3 and i ^ j ^ k. Thus, for the 
projection of spin vectors on the X, Y and Z axis we find 

523^0 = ^ [v{V2 + v;vi + v;v^ + viVz\ , (3.32a) 

^31,0 ^ _:L \y*y^ _ y*y^ + y*y^ _ y*y^ ^ (3.32b) 

Acar 

^'''° " 2^ [^1*^' " ^2*^2 + y^Vz - VIV^ ■ (3.32c) 

The chronometric invariant spin tensor takes the form 

ff^ = {S.,s^S^^^^'\ (3.33) 

and the projection of the spin vector on k axis is defined by 

oo 

Su= I SUf^^dxdydz. (3.34) 



From (3.7) we also write the equations for the invariants 5 — "00, P — iipj^ip and A — ipj^^^ip 



Sa-2gAo = 0, (3.35a) 

Po - 2(m - P) Ao = 0, (3.35b) 

io + 2(m-I?)Po + 2^^5*0 = 0, (3.35c) 
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where Sq = tS, Pq = tP, and Aq = tA, leading to the foUowing relation 

S^ + P'^ + A^ = C^/t^, C^ = const. 



(3.36) 



Let us now solve the Einstein equations. To do it we first write the expressions for the components of the energy- 
momentum tensor explicitly. Using the property of flat space-time Dirac matrices and the explicit form of covariant 
derivative V^, one can easily find 



T° = m5 - F{I, J)+e 

Tl = T^ = T| ^VS + gP- F{I, J) - p. 
Summation of Einstein equations ( 3.22a ), ( 3. 22b ), ( 3. 22c ) and ( |3.22d ) multipHed by 3 gives 



(3.37) 



T 3 
7=2' 



(3.38) 



For the right-hand-side of (3.38) to be a func tion o f t only, the solution to this equation is well-known Bs]. As we 
see in the n ext section, the right-hand-side of ( 3.38| ) is indeed a function of r. Given the explicit form of Ln and Lint 
from ( 3.38| ) one finds the concrete solution for r in quadrature. 

Let us express a, b, c through r. For this we notice that subtraction of Einstein equations (3.22b) and ( [3.22a ) leads 
to the equation 



with the solution 



Analogically, one finds 



at be d f a 
ac be dt 



a ^ f f dt 
-.AexpU/- 



= 0. 



£'2exp X 



dt 



Di — const., Xi — const. 



■ dt 
Dsexp X3 

C \ .1 T 



(3.39) 



(3.40) 



(3.41) 



where D2, -D3, X2, X^ are integration constants. In view of ( |3.23 ) we find the following functional dependence between 
the constants Di, D2, -D3, Xi, X2, X^: 



D2^DiD3, 



X2 — Xi + X3 



Finally, from (3.40) and (3.41) we write a{t),b{t), and c{t) in the explicit form 

a{t) = (i?2^3)'/Vi/3exp 
b{t) = iD^^D3)^/\^/h^v 
c{t) = (i?ii?2)-i/Vi/3exp 



2X1 + X3 r dt 

3 J W) 

X1-X3 f _dt_ 

3 J r{t) 
Xi + 2X3 f dt 

Tit) 



(3.42a) 
(3.42b) 
(3.42c) 



Thus the system of Einstein's equations is completely integrated. 



Defining Hubble constant in analogy with a FRW universe from ( 3.42 ) we obtain 



a," 
a,- 



Y, 



3t 



J = 1,2,3, 



or a generalized one 

H={Hi+H2+H3)/3 = t/3t 
Here ai — a, 0,2 — b, 03 = c. The decelaration parameter given by 



(3.43) 



(3.44) 
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q = 



RR 

~R^ 



for a FRW universe with R being the scale factor can also be generalized for the B-I space-time to obtain 






V a, / V a," 



Inserting ( 3.42 ) into ( 3.46 ) one obtains 



1i = 



T - 2f ^ - Y,T + Y/ 

f2 + 2y,f + y2 



a, / V a, 



i = 1,2,3. 



Let us now go back to the Einstein equation (3.8). Taking the divergence of Einstein equation we obtain 









which in our case reads 



ro + I(T°-Ti) = o. 



Putting T^ and Tl into ( |3.49| ) we obtain 



e + (e +p)- + (m - V)So - GPq = 0, 

T 



(3.45) 



(3.46) 



(3.47) 



(3.48) 



(3.49) 



(3.50) 



where 5*0 = tS and Pq — tP. From (3.35a) and (3.35b) we have (m — I')S'o — GPq ~ 0. Further taking into account 
the equation of state, i.e., p = C^ we lind 



de 



(i + C)e 



with the solutions 



e = 



_£o_ 



dr 

T 



v = 



0, 



Ceo 

t1+C 



(3.51) 



(3.52) 



where eq is the integration constant. Note that the relation (3.52) holds for any combination of the material field 
Lagrangian, e.g., spinor or sca lar or interacting spinor and scalar fields. Thus we see that the right-hand side of ( |3.38 ) 
is a function of t only. Then (3.38|), multiplied by 2f can be written as 



27- f = \i{K{Tl+T^) - 2A)t]t = ^[T)i 



Solution to the equation (3.53) we write in quadrature 

dT 



J^{T)d7 



(3.53) 



(3.54) 



Given the explicit form of F(I, J), from ( |3.54 ) one finds concrete function T(t). Once the value of r is obtained, one 
can get expressions for components "ipjit), j — 1,2,3,4. Thus the initial systems of Einstein and Dirac equations 
have been completely integrated. 

Further we will investigate the existence of singularity (singular point) of the gravitational case, that can be done 
investigating the invariant characteristics of the space-time. In general relativity these invariants are composed from 
the curvature tensor and the metric one. Contrary to the electrodynamics, where there are two invariants only 
( Ji = Ffj^^F'^'^ and J2 = •F^^F^'^), in 4-D Riemann space-time there are 14 independent invariants. They are p6[ 



Ii=R 

h = Rp.iyR'^"^, 

h = Rafipi'R"' ^^ 1 
h = *Raf3fiiyR°' ^^ 1 



(3.55a) 
(3.55b) 
(3.55c) 
(3.55d) 
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h = R" R^^Rafiff}^, 
h = R°' R'^^-^Rafipi', 

h = R°' ^''RafiijpR"''^u^ 

/g ~ -kR"" ^^ Ra(3crpR'^''^,y, 

ha = RaR""^ RlivR'fi, 

111 == RuRp,f°'Raa Rp > 

h2 = R'i*R'"^p^R^£^'' Rp , 

/i3 - R'-^f^iA"^,, + R'^^RPR':,Ry^), 

/i4 = *R^,,rA''^,. 



'■afj 



piV 



(3.55e) 
(3.55f) 
(3.55g) 
(3.55h) 
(3.55i) 
(3.55J) 
(3.55k) 
(3.551) 
(3.55m) 
(3.55n) 



where ^"^^^ = AR^R^R-^R^ + ^R'^Rp,RPRZ and ^i?, 



■p^-fj.^-CT^'-i' """ ■'^"aPp.v — 2 <^P<^P p.i' — 2 '^PP'^ O-P ' ^'aP ~ 2 '^P'^P 

with Eaf)p,i, = ^/—gSa/Siii' and E"^'^'^ = -^==e°'^^'^ . Here Sappu is the totally antisymmetric Levi-Civita tensor with 
£oi23 — 1- Instead of analysing all 14 invariants mentioned above, one can confine this study only in 3, namely the 
scalar curvature /i = R, I2 = R^^j^ii' and the Kretschmann scalar J3 = RafSpuR"^^'^ ■ At any regular space-time point, 
these 3 invariants I — I, I2, I3 should be finite. Let us rewrite these invariants in details. 
For the Bianchi I metric one finds the scalar curvature (see appendix) 



h=R 



f — abc — bca — cab 



(3.56) 



Since the Ricci tensor for the Bianchi I metric is diagonal, the invariant I2 = Rp^R^'^ = RI^Rp is a sum of squares of 
diagonal components of Ricci tensor, that i.e.. 



/2 = 



with 



\Rlf 

Rq ^ - 

Rl = - 
Rl = - 
Rl = - 



{R\) 



[Rlf 



Ri 



abc + abc + abc 



T 




abc + abc + 


abc 


T 




abc + abc + 


abc 


T 




abc + abc + 


abc 



(3.57) 

(3.58a) 
(3.58b) 
(3.58c) 
(3.58d) 



Analogically, for the Kretschmann scalar in this case we have I3 = R^'^ adR'^^pv^ a sum of squared components of all 
nontrivial i?^'' 



Q/3- 



h = A 



R 01 



i?^^n,) +(i?%3) +(i? 12 



4 
72 



{r'\,) 

{dbcf + {abcf + [abcf + [abcf + {abcf + [abc) 



n23 



R 



31 



31 



ibc. 



From (3.42) we have 



a, = ^,ri/3exp[(B,/3) f r-^dt], 

Bt + la, 1 o Q ^ 

3 T 

.. _ jB, + 1){B, ~ 2) a, 
"'- 9 r2' 



(3.59) 

(3.60a) 
(3.60b) 
(3.60c) 
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i.e., the metric functions a,b,c and their derivatives are in functional dependence with r. As we see from (3.60), at 
any space-time point, where r = the invariants /i, I2, I3 become infinity, hence the space-time becomes singular at 
this point. 

4. ANALYSIS OF THE RESULTS 

In this section we shall analyze the general results obtained in the previous section. In the subsections follow we 
will study the system with linear and nonlinear scalar fields respectively. 

A. Linear spinor field in B-I universe 

In this subsection we study the linear spinor field in B-I universe. The reason for getting the solution to the self- 
consistent system of equations for the linear spinor and gravitational fields is the necessity of comparing this solution 
with that for the system of equations for the nonlinear spinor and gravitational fields that permits clarifications of 
the role of nonlinear spinor terms in the evolution of the cosmological model in question. 

In this case we get explicit expressions for the components of spinor field functions and metric functions: 

^i(t) = (Ci/V^)exp [-imt], (4.1a) 

Mt) - (C2/V^)cxp [-imt], (4.1b) 

Mt) - (C3/V^)exp [tmt], (4.1c) 

^|J^{t) = (C4/V^)exp [imt], (4.1d) 



with Ci, C2, C3, Ci being the integration constants. On the other hand from (3.35) we find 

r. Co 



(4.2) 



where Co is an integratio n con stant and related to the previous ones as Cq = Cf + C2 —C^ — C^. For the components 
of the spin current from (|3.27) we find 



f = ^ [C! + ci + ci + ci] , (4.3a) 

f = — [C1C4 + C2C3]cos(2mt), (4.3b) 
ar 

f = ^ [C1C4 - C2C3]sm{2mt), (4.3c) 

OT 

f = — [C1C3 - C2C4]cos{2mt), (4.3d) 



CT 

whereas, for the projection of spin vectors on the X, Y and Z axis we find 

S^^^°^^[CiC2 + C3C4], (4.4a) 

S^^-° = 0, (4.4b) 

S'''' = ^[C!~Ci + Ci-Cl]. (4.4c) 

From ( 3.29| ) we find the charge of the system in a volume V 

Q = [Cf + Ci + Ci + Cj] V. (4.5) 

Thus we see, the total charge of the system ina finite volume is always finite. 

Let us now determine the function r. In absence of perfect fluid for the linear spinor flcld we have 

To° ^ mS, Tl = r| = T| = 0. (4.6) 

Taking ( |4.6D into account, for r we write 

f = M - 3Ar (4.7) 
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with the solutions 



(1/3A) [M - (7isinh(V^3Ai)], 
(l/3A)[M-g2sin(V3Ai)], 



A<0, 
A = 0, 
A>0. 



(4.8) 



where M — ^kuiCq and yi 



„,u,,„^y, i^^i^ yj^, j/Oj 9ii 92 are the constants. Let us now analyze the solutions obtained. 
First we study the case when A = 0. It can be shown that BQ] 



yl - 2My„ = (X^ 
This means that the quadratic polynomial (l/2)Afi^ - 



turns into zero at i = 
in this case we have 



(4.9) 

T(t) in case of A = 
ti_2 = —yi/M ± ^J{yi/MY — 2yo/M and the solution obtained is the singular one. At t — > cx) 



~a:iX3 + a:|)/3>o. 

Vii + yo = possesses real roots, i.e. 



r(t) K, —KmCot 



ait) w b{t) w c{t) w t 



2/3 



which leads to the conclusion about the asymptotical isotropization of the expansion process for the initially anisotropic 
B-I space. Thus the solution to the self-consistent system of equations for the linear spinor and gravitational fields is 
the singular one at the initial time. In the initial state of evolution of the field system the expansion process of space 
is anisotropic, but at t — > oo the isotropization of the expansion process takes place. As one sees the components 
of spin current and projections of spin vector are singular at space-time points ii_2 where r vanishes. A qualitative 
picture of this case has been given in Fig. 1. 

For A < we see that the solution is singular a,t t = to = (1/v— 3A)arcsinh(A//gi) and the isotropization of the 
expansion process takes place as i — > oo. Note that the izotropization process in this case is rather rapid [cf. Fig 2]. 

For A > we have the oscillatory solutions [cf. Fig. 3]. Taking into account that r is a non- negative quantity, 
it can be shown that the model has singular solutions at i = (4fc -t- l)7r/2'\/3A, fc = 0, 1, 2, 3, .... with M = q2. For 
M > q2 we have t that is always positive definite, i.e., the solutions obtained are regular at each space-time point. 
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FIG. 1. Perpective view of r for linear spinor field in absence of A term. 
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FIG. 2. Perpective view of r for linear spinor field with A < 0. 
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FIG. 3. Perpective view of r for linear spinor field with A > 0. 



B. Nonlinear spinor field 

Let us now go back to nonlinear case. We consider the following forms of nonlinear term: I. Ln = F{I): II. 
Ln = F{J); III. Ln = F{K±) with K± = I ± J. 

I. Let us consider the case when Xn = F{I)- From ( 3.35| ) we find in this case we find 



S = — , Co — const. 



(4.10) 



Note that in this case we denote the constants in the same way as we did it for linear case, but the constants in these 
cases are not necessarily identical. Spinor field equations in this case read 



Vi+i{m~V)Vi =0, 
V2+iim~V)V2 = 0, 
V3-i{m~V)V3^0, 
Vi - i{m - V)Vi = 0. 



(4.11a) 
(4.11b) 
(4.11c) 
(l.lld) 



As in the considered case when Ln = F depends only on S, from ( 4.1C| ) it follows that F{I) and V are functions of 
r only. Taking this fact into account, we get explicit expressions for the components of spinor field functions 

(4.12a) 

(4.12b) 

(4.12c) 

(4.12d) 

with Ci, 6*2, C3, C4 being the integr ation constants and are related to Co as Co = Cf + C2 ~ C^ — C^- For the 
components of the spin current from (3.27) we find 



^1 {t) = (Ci / V7)exp [-ij{m- V)dt] , 
Mt) = (C2/V7)exp[-z /"(m-I?)di], 
i^z{t) = (C3/\/7)exp [i j{m- V)dt], 
Mt) = (C4/V^)exp [i [m - V)dt], 



f^-[cf + cl + cl + cl], 

j^ = — [C1C4 + C2C3]cos[2 /(m - V)dt], 
at'- 'J 

f ^^[CiCi - C2C3]sm[2 J {m - V)dt], 
f = — [C1C3 - C2C4]cos[2 f{m - V)dt], 



(4.13a) 
(4.13b) 
(4.13c) 
(4.13d) 
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whereas, for the projection of spin vectors on the X, Y and Z axis we find 



c23,0 
^31,0 ^ 
^12,0 ^ 



^[CiC2 + C3C4], 



2abT 



[cl~cl + cl-cl]. 



(4.14a) 
(4.14b) 
(4.14c) 



We now study the equation for r in details choosing the nonhnear spinor term as F{I) — A/^"/^-' — AS*" with A 
being the couphng constant and n > 1. In this case for r one gets 



T 



r = (3/2)K[mCo + \{n - 2)C^7 
The first integral of the foregoing equation takes form 

f2 = 3K[mCoT - ACo"/r"-2 + g^ 



3A7 



3Ar^ 



(4.15) 



(4.16) 



Here g^ is the integration constant that is positive defiend and connected with the constants Xi as g^ = {Xi+XiX^ + 
X|)/9kCo p3]. The sign Cq is determined by the positivity of the energy-density Tg of linear spinor field, i.e., 

To° = mCo/r > 0. (4.17) 



It is obvious from ( 4.17 ) that Cq > 0. Now one can write the solution to the equation ( 4.16 ) in quadratures: 



/ 



yKfrnCoT" 



„2^n-2 



XCll 



At" 



= %/3t 



(4.18) 



The constant of integration in (4. IS) has been taken to be zero, as it only gives the shift of the initial time. Let us 
study the properti es of solution obtained for different choice of n, A and A. Firts we study the case with A = 0. 
For n > 2 from (fTq) one gets 



-(t) It^oo- (3/4)KmCoi2. 



(4.19) 



It leads to the conclusion about isotropization of the expansion process of the B-I space-time. It should be remarked 
that t he iso tropization takes place if and only if the spinor field equation contains the massive term [cf. the parameter 
m in (4.18|)]. This is not the case for a massless spinor field, since from (4.18) we get 



T(t) 



\/'inCag^ t. 



(4.20) 



Substituting (4.20) into (3.42) one comes to the conc lusio n that the functions a{t),b{t), and c(t) are different. 
Let us consider the properties of solutions to Eq. ( [4.15 ) when i — > 0. For A < from ( 4.18| ) we get 

T{t) = [{3/A)n\\X\C'o'] '^"t'Z" ^ 0, (4.21) 

i.e. solutions are singular. For A > 0, from ( [4.1^ ) it follo ws that r = cannot be reached for any value of t as in 
this case when the denominator of the integrand in ( [l.lSl) becomes imaginary. It means that for A > there exist 
regular solutions to the previous system of equations p9|. The absence of the initial singularity in the considered 
cosmological solution appears to be consistent with the violation for A > of the dominant energy condition in the 
Hawking-Penrose theorem [Esl which reads as follows: 

THEOREM. A space-time M. cannot he causally, geodesically complete if the GTR equations hold and if the 
following conditions are fulfilled: 

1. The space-time Ai does not contain closed time- like lines. 

2. The conditions (dominant energy condition) 



Too + 2^11 + 2122 + Tas > 0, 
Too + Tn > 0, 
Too + T22 > 0, 
Too + T33 > 0, 



(4.22a) 
(4.22b) 
(4.22c) 
(4.22d) 



on the equations of state are fulfilled, where Too is the energy density and Tn, r22, and T^s are three principal values 
of pressure tensor. 
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3. On each time-like or null geodesic, there is at least one point for which 



(4.23) 



where Ka is the tangent to the curve at the given point and where the brackets on the subscripts imply antisym- 
metrization. 

4. The space-time M contains either (a) a point P such that all diverging rays from this point begin to converge 
if one traces them back into the past, or (6) a compact space-like hypersurface. 

To prove that in the case considered the dominant energy condition violates, we rewrite (4.22) in the following 
form: 

(4.24a) 
(4.24b) 
(4.24c) 
(4.24d) 



rpQ 
^0 


> T^a^ 


+ T^b' 


+ Tic' 


^0 


> Tla" 






J-0 


> nh\ 








> Tic'. 







Let us go back to the energy density of spinor field. From 






mCo XC^ 



(4.25) 



follows that at 



^„_i ^ ACq 



n — 1 



(4.26) 



the energy density of spinor field becomes negative. On the other hand we have 



T^ = Ti = Tl = 



\{n - l)Co" 



>0 



(4.27) 



for any non-negative value of r. Thus, we see all four conditions in (4.24) violate, i.e., the absence of initial singularity 
in the considered cosmological solution appears to be consistent with the violation of the dominant energy condition 
in the Hawking-Penrose theorem. 

Let us consider the Heisenberg-Ivanenko equation ^^ setting n = 2 in ( 4.15 ). In this case the equation for T(i) 
does not contain the nonlinear term and its solution coincides with that of the linear one. The spinor field functions 
in this case are written as follows: 



y _ _^ -imt j^iiXCo / B 



V2 



JT 
C2 _ 



imt r^4iXCo / B 



Q 



3_ imt y-iiXCo/B 



Ci 



Y _ J^ imty-Ai\Co/B 



(4.28a) 
(4.28b) 
(4.28c) 
(4.28d) 



where Z = L^^j , B = M{ti — ^2), and ti^2 — ~yi/M ± •^/(j/i/Af)^ — 2yQ/M are the roots of the quadratic equation 

Mt^ -|- 2yit + 2j/o = 0. As in the linear case, the obtained solution is singular at initial time and asymptotically 
isotropic as i — > 00. 

We no w stu dy the properties of solutions to Eq. (4.15) for 1 < n < 2. In this case it is convenient to present the 
solution (4.18) in the form 



dr 



mT-AT2-"C^*-i+c,2 



= v 3kCo t 



(4.29) 



As t ^ 00, from ( |4.29| ) we get the equality ( 4.19| ), leading to the isotropization of the expansion process. If m = 
and A > 0, T(i) lies on the interval 



< T{t) < (5VAQ-I) 



-ni/Cz-T- 
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If m = and A < 0, the relation ( 4.29| ) at i — > oo leads to the equality 



(4.30) 



Substituting ( 4.30 ) into ( 3.42 ) and taking into account that at t 

dt n(3K|A|n2Co")i/" 



[n - 2)22/" 



i-2/n+l ^ 



due to —2/n + 1 < 0, we obtain 



a{t) ~ h{t) - c(i) - [T(i)]'/^ - t'/^" 



(4.31) 



This means that the solution obtained tends to the isotropic one. In this case the isotropization is provided not by 
the massive parameter, but by the degree n in the term Lj^ = A5". Equation (4.29) implies 



'it) 



\/3KCog^ t ^ 0, 



(4.32) 



which means the solution obtained is initially singular. Thus for 1 < n < 2 there exist only singular solutions at 
initial time. At i ^ cxd the isotropization of the expansion process of the B-I space takes place both for m ^ Q and 
for m = 0. 



Finally, let us stu dy th e properties of the solution to the equation ( 4.15 ) for < n < 1. In this case we use the 
solution in the form ( 4.2£ ). Since now 2 — n > 1, then with the increasing of T(i) in the denominator of the integrand 
in (4.29) the second term At2~"Cq ~^ increases faster than the first one. Therefore the solution describing the space 
expansion can be possible only for A < 0. In this case at i — *■ oo, for m = as well as for m ^ 0, one can get the 



asymptotic representation (4.3C) of the solution. This solution, as for the choice 1 < n < 2, provides asymptotically 
isotropic expansion of the 15-1 space-time. For t ^ in this case we shall get only the singular solution of the form 

(HI)- 

For a nonzero A term we study the following situations depending on the sign of A and A. 
case 1. A = — e^ < 0, A > 0. In this case for n > 2 and t ^ oo we find 

r(i) « 6^"=* (4.33) 

Thus we see that the asymptotic behavior of r does not depend on n and defined by A - term. From (|3.42) it is 



obvious that the asymptotic isotropization takes place. 

From ( 4.1^ ) it also follows that r cannot be zero at any moment, since the intigrant turns out to be imaginary 
as T approaches to zero. Thus the solution obtained is a nonsingular one thanks to the nonlinear term in the Dirac 
equation and asymptotically isotropic. As it has been noted earlier, the absence of initial singularity in the considered 
cosmological model results the violation of the dominant energy condition. 

case 2. A > and A > 0. For n > 2 ( 4.1g ) admits only nonsingular oscillating s oluti ons, since r > and bound 
from above. Consider the case with n = 4 and for simplicity set m — Q. Then from ( 4.18 ) one gets 



T[t) 



1 



kCqto 



K CqTq 



4AAC^ sin2V3Ai 



1/2 



(4.34) 



For a massive spinor field with A > and A > and n = 10 a perspective view of r is shown in FIG. 4. The period 
for the massive field is greater than that for the massless one. As it occurs, the order of nonlinearity (n) has a direct 
effect on the period (the more in n the less is the period) 



1.55 
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0.95 
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FIG. 4. Perspective view of r showing the initially nonsingular and oscillating behavior of the solutions. The continuous and 
dash lines correspond to the massive and massless spinor field, respectively. 

case 3. A < and A < 0. The solution is singular at initial moment, that is 



[^-3An2Co"/4i]2/" (4.35) 

and at i — > oo asymptotic isotropization takes place since 



lim r : 



lim T« e^-^"^*. (4.36) 

t—*oo 



case 4. A > and A < 0. Solution is initially singular as 



limr«[^-3An2Co"/4i]'/" (4.37) 

and bound from the above, i.e., oscillating, since 

lim T « sin\/3At. (4.38) 

II. We study the system when L^ = F{J), which means in the case considered V — 0. Let us note that, in the 
unified nonlinear spinor theory of Heisenberg, the massive term remains absent, and according to Heisenberg, the 
particle mass should be obtained as a result of quantization of spinor prematter p9| ] . In the nonlinear generalization 
of classical field equations, the massive term does not possess the significance that it possesses in the linear one, as it 
by no means defines total energy (or mass) of the nonline ar fie ld system. Thus without losing the generality we can 
consider massless spinor field putting m = 0. Then from ( 3.35| ) one gets 

P(i) = — , Do = const. (4.39) 

r 

The system of spinor field equations in this case reads 

Vi - gv3 = 0, 

V2 - QVi = 0, 

1/3 + gvi = 0, 
1/4 + gv2 = 0. 

Defining U{a) — V{t)^ where u = jQdt, we rewrite ( 4.40| ) as 

U[-U3 = 0, 
t/^ - C/4 = 0, 
t/^ + C/i = 0, 

K + U2 = 0, 

where prime (') denote differentiation with respect to a. Differentiating the first equation of system ( 4.41 ) and taking 
into account the third one we get 



(4.40a) 


(4.40b) 


(4.40c) 


(4.40d) 


(4.41a) 


(4.41b) 


(4.41c) 


(4.41d) 



which leads to the solution 



Analogically for U2 and f/4 one gets 



U[' + Ui= 0, (4.42) 



Ui = Die"" + iD^e- 
U3 = iDie"' + D-se- 



U2 = D2e"' + lD^e^ 
Ui = iL'se*'" + Die' 



where Di are the constants of integration. Finally, we can write 
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i^i = ^(Die" 



iDse' 



Putting ( 4.43| ) into the expressions ( 4.3£ ) one comes to 

Dn = 2[Dl + Dl-Dl- DJ). 



For the components of the spin current from (3.27) we find 

,f = -[Dl + Dl + Dl + Dl], 



/ ^—[D2D3 + DiDi\cos[2 [gidt], 
f == ^[^2^3 - DiD4]sm[2 fgidt], 
f ^^[DiDs - D2Di\cos[2 I gi)dt], 



whereas, for the projection of spin vectors on the X, Y and Z axis we find 

2 



S 



23,0 



^^^[DiD2 + D3Di], 

^31." = 0, 
^12,0 ^ 1 



2abT 



[Dl-Dl+Dl-Dl]. 



Let us now estimate t using the equation 



f/r = 3KX{n- l)P 



2n 



where we chose Ln = XP . Putting the value of P into (4.46) and integrating one gets 



f2 = -3/cAZ?2V-2"+y2, 



(4.43a) 
(4.43b) 
(4.43c) 
(4.43d) 



(4.44a) 
(4.44b) 
(4.44c) 
(4.44d) 

(4.45a) 
(4.45b) 
(4.45c) 

(4.46) 
(4.47) 



where y^ is the integration constant having the form y^ — {X^ + XiX^ + X^)/3 > 0. The solution to the equation 
( [4.47 ) in quadrature reads 



dr 



^-3kA7:>^"t2-2" + y2 



(4.48) 



Let us now analyze the solution obtained here. As one can see the case n = 1 is the linear one. In case of A < for 
n > 1 i.e. 2 — 2n < 0, we get 



it) |,^o« [{V^^X\D^n)t]'/", 



and 



r t- 



V3Kj/^i. 



This means that for the term Ln considered with A < and n > 1, the solution is initially singular and the space-time 
is anisotropic at t -^- cxd. Let us now study it for n < 1. In this case we obtain 



|t^o~ \piK\p-t 
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and 

The solution is initially singular as in the previous case, but as far as l/n > 1, it provides an asymptotically isotropic 
expansion of B-I space-time. The analysis for A 7^ completely coincides with those for F = A5" with m = 0. 

III. In this case we study Ljv — F{I, J). Choosing 

Ln = FiK±), K+ = I + J = I, = -Ia, K^=I-J = It, (4.49) 

in the case of massless NLSF we find 

V = 2SFk^, g = ±2PFk±, Fk^ = dF/dK±. 

Putting them into ( 3.35| ) we find 

Sl±P^^D±. (4.50) 



Choosing F — \K!j_ from ( 3.3q ) we get 

r = 3KX{n - 1) Dl r^-^", (4.51) 

with the solution 

= i, (4.52) 



where g^ = (X^ + X1X3 + X|)/3. Let us study the case with A < 0. For n < 1 from (3.33) one gets 

T{t)\t^o~9t~^0, (4.53) 

i.e., the solutions are initially singular, and 



rit) |,^oo« [^(3«|A|i?S)t]i/", (4.54) 

which means that the anisotropy disappears as the Universe expands. In the case of n > 1 we get 

T(i) \t^o~ t'/" - 0, 
and 

T{t) It^oo- gt, 

i.e., the solutions are initially singular and the metric functions a{t),b{t), and c{t) are different at i ^ cxd, i.e., the 
isotropization process remains absent. For A > we get that the solutions are initially regular, but it violates the 
dominant energy condition in the Hawking-Penrose theorem ES] . Note that one comes to the analogical conclusion 
choosing Lw = XS^'^P^''. 

C. Analysis of the results obtained when the B-I Universe is filled with perfect fluid 

Let us n ow analyze the system filled with perfect fluid. In absence of other matter, i.e., spinor field, in this case 
from ( [3.3§| ) we find 

i|ii-f££, (4.55) 

with the first integral 



V3Keor(i-C)+C, (4.56) 



where C is an integration constant. From (4.56) one estimates 
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T(xt'^, for C = 0, (dust), 
Tcxt^/^, for ^=1/3, (radiation), 
T ex t^'^, for C = 2/3, (hard universe), 
T oc t, for C = 1> (stiff matter). 

A perspective view of these solutions are given in Fig. 5. 



(4.57a) 
(4.57b) 
(4.57c) 
(4.57d) 



io 




io 



FIG. 5. Perpective view of r when B-I universe is filled with perfect fluid only. The lines from left to right at the upper 
corner correspond to dust ((" = 0), radiation {( — 1/3), hard universe (^ = 2/3) and stiff matter {( — 1), respectively. 



Let us now consider the system as a whole with the nonlinear term being L^ — XS" . f n this case we get 

dT 



I 



y/mCoT - ACJ/t("-2) + eor(i-«) + g^ 



±V3Kt. 



(4.58) 



As one can see in the case of dust (^ — 0) the fluid term can be combined with the m assive one, whereas in the case 
of stiff matter (^ — 1) it mixes up with the constant. Analyzing the equation ( 4.58) one concludes that in presence 
of spinor field perfect fluid plays a secondary role in the evolution of B-I universe. 



5. CONCLUSION 

Within the framework of the simplest nonlinear model of spinor field it has been shown that the A term plays very 
important role in Bianchi-I cosmology. In particular, it invokes oscillations in the model which is not the case when 
A term remain absent. It should be noted that regularity of the solutions obtained by virtue of A term, specially for 
the linear spinor field does not violate dominant energy condition, while this is not the case when regular solutions 
are attained by means of nonlinear term. Growing interest in studying the role A term by present day physicists of 
various discipline witnesses its fundamental value. For details on time depending A term one may consult |3q| and 
references therein. 



[1 
[2 

[s; 

[4 

[6; 
[7; 

[9 

[to; 
[11 



A.F. Ranada, in A.O. Barut (ed.) "Quantum theory, groups, fields and particles", 271 

D. Ivanenko, Soviet Physics, 13, 141-150 (1938) 

D. Ivanenko, in Einstein Centenay Jubilee Volume. 

V. Rodichev, Soviet Physics JETP, 13, 1029-1033 (1961) 

H. Weyl, Physical Review, 77, 699-701 (1950) 

R. Utiyama, Physical Review, 101, 1596-1607 (1956) 

T.W.B. Kibble, Journal of Mathematical Physics, 2, 212-221 (1961) 

D.W. Sciama, Festschrift for Infeld, Pergamon Press, 1960, 415-439. 

A.F. Ranada and M. Soler, Journal of Mathematical Physics, 13, 671-676 (1972). 

F.W. Hehl, P. von der Heyde and G.D. Kerlick, Review of Modern Physics, 43, 393-416 (1976) 

W. Heisenberg, Physica 19, 897-908 (1953) 



291. Reidel (1983) 



Spinor field in Bianchi type-I universe: regular solutions 21 

[12] W. Heisenberg, Review of Modem Physics, 29, 269-278 (1957)]. 

[13] D.J. Gross and A. Neveu, Pliysical Review D, 10, 3235-3253 (1974) 

[14] L. Parker, Pliysical Review 183, 1057 (1969). 

[15] L. Parker, Physical Review D 3, 346 (1971). 

[16] C.W. Misner, Astrophysics Journal 151, 431 (1968). 

[17] B.L. Hu and L. Parker, Physical Review D 17, 933 (1978). 

[18] G.F.R. Ellis, Journal of Mathematical Physics, 8, 1171 (1967) 

[19] J.M. Stewart and G.F.R. Ellis, Journal of Mathematical Physics, 9, 1072 (1968) 

[20] V.A. Belinskii, I.M. Khalatnikov and E.M. Lifshitz, Advances in Physics 19, 525 (1970). 

[21] K.C. Jacobs, Astrophysics Journal 153, 661 (1968). 

[22] L.P. Chimento and M.S. MoUerach, Physics Letters A 121, 7 (1987). 

[23] M.A. Castagnino, CD. El Hasi, F.D. MazziteUi and J. P. Paz, Physics Letters A 128, 25 (1988). 

[24] B.L. Hu, Physical Review D 18, 968 (1978). 

[25] P.G. Miedema and W.A. van Leeuwen, Physical Review D 47, 3151 (1993). 

[26] H.T. Cho and A.D. Speliotopoulos, Physical Review D 52, 5445 (1995). 

[27] G.N. Shikin, Preprint IPBRAE, Academy of Science of USSR 19, 21 p. (1991). 

[28] Yu.P. Rybakov, B. Saha and G.N. Shikin, PFU Reports, Physics, 2, (2), 61 (1994). 

[29] Yu.P. Rybakov, B. Saha and G.N. Shikin, Communications in Theoretical Physics 3, 199 (1994). 

[30] B. Saha and G.N. Shikin, Journal Mathematical Physics 38, 5305 (1997). 

[31] R. Alvarado, Yu.P. Rybakov, B. Sah a and G.N. Sh ikin, JINR Preprint E2-95-16, 11 p. (1995), Communications in 



Theoretical Physics 4, (2), 247 (1995), ;r-qc/9603035 



[32] R. Alvarado, Yu.P. Rybakov, B. Saha and G.N. Shikin, Izvestia VUZob 38,(7) 53 (1995). 

[33] B. Saha, and G.N. Shikin, General Relativity and Gravitation 29, 1099 (1997). 

[34] Bijan Saha, and G.N. Shikin, On the role of A-ter m in the evolut ion of Bianchi-I cosmological model with nonlinear spinor 



field, (to be published in PFU Reports. Physics.) gr-qc/010205£ 



[35] Bijan Saha, Dirac Spinor in Bianchi-I Universe with t ime dependent Gravitational and Cosmological Constants, (to be 



pubhshed in Modern Physics Letters A) gr-qc/0009002 



[36] A. Einstein, Sitz. Ber. Preuss. Acad. Wiss. (1917) and (1919) [English translation; The Principle of Relativity (Methuen, 

1923, reprinted by Dover, New York, 1924), pp. 177 and 191. 
[37] L. Bianchi, Lezioni Sulla Teoria dei Gru ppi Continui F initi Transjormazioni (Spoerri, Piza, 1918). 



[38] M. Tsamparhs and P.S. Apostolopoulos, ;r-qc/0006083 



[39] Silvan S. Schweber, An Introduction to Relativistic Quuantum Field Theory (Row, Peterson and Co., NY, 1961). 

[40] V.B. Berestetski, E.M. Lifshitz and L.P. Pitaevski, Quantum Electrodynamics (Nauka, Moscow, 1989). 

[41] V.A. Zhelnorovich, Spinor theory and its application in physics and mechanics (Nauka, Moscow, 1982). 

[42] D. Brill, and J. Wheeler, Review of Modern Physics 29, 465 (1957). 

[43] Ya.B. Zeldovich, and LD. Novikov, Structure and evolution of the Universe (Nauka, Moscow, 1975). 

[44] N.N. Bogoliubov, and D.V. Shirkov, Introduction to the theory of quantized fields (Nauka, Moscow, 1976). 

[45] E. Kamke, Differentialgleichungen losungsmethoden und losungen (Leipzig, 1957). 

[46] N.V. Mitskevich, A. P. Efremov and A.I. Nesterov, Dynamics of fields in general relativity (Energoatomizdat, Moscow, 

1985). 
[47] D.D. Ivanenko, Introduction to Nonlinear Quantum Field Theory (Foreign Literature, Moscow, 1959). 
[48] S.W. Hawking, and R. Penrose, Proceedings of Royal Society, London, A 314, 528 (1970). 
[49] W. Heisenberg, Introduction to the unified field theory of elementary particles (Interscience Publ., London. 1966). 



